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RESUMEN 

Cuatro transitos del plancta extrasolar HAT-P-23b fucron observados recien- 
temente con el telescopio de 0.36 m del Observatorio de la Universidad de 
Monterrey. Las cuatro curvas de luz fueron exitosamente combinadas para 
obtcner una rcsultante de mayor calidad. Esta fue modclada utilizando un 
metodo de Monte Carlo para obtener los parametros esenciales que carac- 
tcrizan al sistema. Asumiendo parametros orbitales como la execntricidad 
e y la longitud del periastro uj rcportados en el articulo del descubrimiento, 
encontramos valores de Rp/R+ = 0.1105+q'ooi3 para la razon del radio del 
planeta al radio estclar, a/R+ = 4.23^^2 para el scmicje mayor normalizado, 
y una inclination orbital de i = 87.9° j^a- Obtenemos un periodo orbital de 
1.2128868 ±0.0000004 dfas (T Q = 2, 454, 852.26542 ± 0.00018 BJD.TDB). 

ABSTRACT 

Four transits of the exoplanet HAT-P-23b were recently observed with the 
0.36 m telescope at the Universidad de Monterrey Observatory. The four 
light curves were successfully combined to obtain a resulting one with re- 
duced scattering per bin. This one was modeled using a Monte Carlo method 
to obtain the essential parameters that characterize the system. Assuming 
orbital parameters such as eccentricity e and longitude of periastron oj from 
the discovery paper, we found values of Rp/R± = 0.1105lo!ooi3 ^ or * ne P iane t- 
to-star radius ratio, a/R± = 4.23^q'j2 for the scaled semimajor axis, and an 
orbital inclination of the system of i = 87.9° ^ \- We also derive an improved 
orbital period of 1.2128868 ±0.0000004 days (T Q = 2, 454, 852.26542 ±0.00018 
BJD.TDB) for the system. 
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1. INTRODUCTION 

HAT-P-23b is a relatively massive (~ 2Mjtjp) transiting exoplanet or- 
biting the GO dwarf star GSC 1632-01396 (m v ~ 12.4) on a close circular 



orbit with a period of 1.212884 ± 0.000002 days (jBakos et al.ll2011l) . In the 

discovery paper, the authors conclude that this exoplanet is an inflated hot 
Jupiter, which has one of the shortest characteristic in-fall times (~ 7.5^ 2 'g 
Myr) before being engulfed by the star. The Rossiter-McLaughlin effect for 
HAT-P-23b was expected to be significant because of the moderately high 
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rotational velocit y of the star ( 8 .1 ± .5 km s -1 ) and transit depth (~ 17 
mmag). However. iMoutou et al.l (j 2 1 lh measured a similar wsini^ (7.8 ± 1.6 
km s _1 ), but estimated a prograde aligned orbit with a projected angle be- 
tween the orbital plane and the stellar equatorial plane of only A = +15° ±22° 
and suggested that the uncertainty in the inclination of the system has some 
impact on the relatively large error bars. Thus it is of importance to at- 
tempt to improve this parameter from light curve observations. It is also 
important to verify the radius of the plan et since the current reported radius 
R p = 1.368 ± 0.090i?j (IBakos et al.ll201ll) cannot be reproduced by the the- 
oretical models developed by iFortnev et al. I d2008h . which suggest a smaller 
size for a planet of this mass. For planets that are larger than expected, mod- 
els generally co nsider an internal sou rce of heating to explain the increase in 
planetary radii ( Laughlin et alJ[2011 ) . Possible inflation mech anisms are tidal 
dissipation in orbit circularization (jBodenheimer et al.ll2001l) . kinetic energy 
transpor t from wind-driven activit y in the atmosphere to dissipation at the 
interior (jGuillot fc Showman! 120021) . and internal heating by Ohmic dissipa- 
tion of currents driven to the interior, produced by magnetohydrodynami c 
mechanisms on the surface ( Batvgin fc Stevenson 2010t Laughlin et al. 2011 ). 

In this paper, we analyze recent transits of the extrasolar planet HAT-P- 
23b in an attempt to improve the essential model parameters that characterize 
this system. We describe a method for reducing and combining the observa- 
tions of each transit to obtain an average light curve from which the transit 
parameters can be obtained. In § El we discuss the methodology for the pho- 
tometric observations and data reduction of the transits. In § [3J we describe 
our method for combining several transit curves in order to obtain a curve 
with reduced noise per bin, and we derive an improved orbital period for the 
system. In § [4j we obtain the transit parameters from a best-fit model. In §[5j 
we compare our results with previous observations and in § [6] we summarize 
our findings. 



2. OBSERVATIONS AND DATA REDUCTION 

We observed four transits of HAT-P-23b using the Universidad de Monter- 
rey (UDEM) Observatory telescope on UT dates 2011 June 3, and on August 
4, 16 and 21. This is a small private college observatory having Minor Planet 
Center Code 720 and is located at an altitude of 689 m, in the suburbs of Mon- 
terrey, Mexico. The data were acquired using a standard Rc-band filter (630 
nm) on a 0.36 m reflector with a 1280 x 1024 pixel CCD camera at 1.0 arcsec 
pixel -1 scale, resulting in a field-of-view of ~ 21.3 x 17.1'. The observations 
were slightly defocused. This technique spreads light over more pixels, which 
allows longer exposure times without saturation and reduces systematic errors 
of focusing light on a few pixels, thus increasing photometric precision. Also, 
on-axis guiding was used to maintain pointing stability. The exposure times 
were 60 s and the images were binned 2 x 2 to facilitate rapid readout (~ 3 s). 
Each observing session lasted at least four hours in order to accommodate the 
transit event and also to cover about one hour before ingress and one hour 
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after egress. Weather conditions were similar during all the observing ses- 
sions (clear and calm), with the August 4th night being a bit more windy and 
turbulent. In order to obtain accurate time stamps for our images, the data- 
taking computer was synchronized with Coordinated Universal Time (UTC) 
at the start of each observing session through the Internet. At the end of 
the observing session, time was independently verified, to the nearest second, 
using WWV radio time signals. 

Standard dark current subtraction and twilight sky flat-field division pro- 
cess were performed on each image for calibration. Aperture photometry was 
carried out on the target star and five comparison stars of similar (± ~ 1.5) 
magnitude. The apertures used varied for each date due to the defocus and 
weather conditions, but were optimized to minimize the scatter of the resulting 
light curves. We found that slightly smaller scatter in the final light curves 
was obtained by averaging the magnitude differences of HAT-P-23 to each 
comparison star individually. This produced smaller scatter than the method 
of ratioing the target star flux to the sum of the fluxes of all comparison stars. 
We found that both methods were essentially equivalent when using a limited 
number of comparison stars with similar magnitudes. Our chosen method 
has the advantage of easily identifying and correcting unsuitable comparison 
star measurements (due to saturated pixel effects, faintness, cosmit rays, etc.) 
and also provides an independent measure of the uncertainty for each point, 
which was similar to the overall data scatter. We estimated the formal error 
for each HAT-P-23 photometric point as the standard deviation of the mag- 
nitude differences to the individual comparison stars, divided by the square 
root of their number (error of the mean). 

After normalizing the target star to the comparison stars and averaging, 
some gradual variations as a function of time were found. This is perhaps 
caused by differential extinction between the transit and comparison stars, 
which generally have different and, in some cases, unknown spectral types. 
Bluer comparison stars are more affected by this atmospheric scattering effect 
than redder stars. Consequently, this variation was removed by using a linear 
air mass-dependent function of the form: 

5m = c(l - X) + b (I) 

where 8m is the magnitude change applied to an individual measurement at air 
mass A, and c and b are best-fit constants necessary to remove the systematic 
effect and obtain a flat (horizontal) line for the out-of-transit baseline portions 
of the light curve. The best-fit model in equation ([TJ is determined from the 
out-of-transit data. All transits have at least one hour's worth of baseline 
observations before the transit ingress and after the transit egress for this 
purpose. This could also have been corrected by introducing an intrinsic color 
term for the stars (using the 2MASS catalogue for example), but our chosen 
empirical procedure yields similar good results and is simpler to implement. 

The resulting light curves for the four observing dates are presented in 
Fig. HI 
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3. LIGHT CURVE ANALYSIS 
3.1. Light Curve Combination 

Transit photometry may present significant scattering due to various noise 
sources, specially in ground observations with small telescopes. Therefore, 
we propose to combine all the individual transit curves to produce a final 
average light curve with reduced noise per bin. This is warranted because all 
four light curves are from independent events, have similar scatter, the same 
exposure times, same number of comparison stars, and are devoid of star-spot 
activity. Combining all four transit light curves into one, should theoretically 
decrease the scatter by a factor of ~ y/n. For our observations, this means a 
decrease by a factor of « 2. This procedure also has the tendency to cancel 
out any systematic effects that may be present in the individual light curves. 
From a computational point of view, this scheme facilitates the analysis of 
a relatively large number of subsequent follow-up photometry data, as the 
parameter space would increase by another mid-transit T c , making the model 
analysis more costly. 

In order to combine the light curves, we first need to co-register them 
by finding the mid-transit time for each event. In order to fit the observed 
transit light curves we first created initial standard model light curves. These 
were constructed n umerically as a t ile-the-star procedure using the Binary 



Maker II software ( Bradstreetl 120051 ) . The initial system parameters used 



were from iBakos et al.1 (j 2 1 lh and the linear limb-darkening function coeffi 



cient (u = 0.610 for our bandpass) fo r this system was tak en from 



(2000). Following a method similar to ISada et al.l (|2012l ) and iTodorov et al 



Claret 



2012J), small adjustments to the duration and depth of the model transits 
were necessary to optimize the fits and extract the best mid-transit time pos- 
sible. This was done by applying small (a few percent) multiplicative factors 
to both the depth and duration of the model transit. Best-fit models were 
obtained by minimizing the x 2 of the data. We used the scatter of the pho- 
tometry to estimate the uncertainties because the scatter is larger than the 
formal errors suggest. 

Once the initial mid-transit times for each system were obtained, a com- 
bined light curve was created by placing all the photometry points in a com- 
mon time reference frame, centered on the mid-transit time for each event, 
and by averaging the data points falling within predetermined bins. This 
method generates a light curve with reduced noise per bin, equally spaced 
data points, and maintains the general shape of the transit curve. However, 
care must be taken in selecting the bin size in order to simultaneously have 
sufficient amount of data in each bin, and also sufficient data points to define 
the critical ingress and egress portions of the light curve. After testing differ- 
ent bin sizes, it was determined that 2-minute bins satisfied both criteria and 
this bin size was adopted for averaging the individual transits HAT-P-23b. 

Further refinements were performed to the individual light curve mid- 
transit times by subjecting the combined light curve to the same x 2 fit to 
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the preliminary model and obtaining general multiplicative factors to both 
the depth and duration of the model transit. These factors were then used 
on the individual light curves to derive improved mid-transit times and gen- 
erate a new combined light curve. The process was iterated until there was 
convergence between the individual mid-transit times and the best-fit multi- 
plicative factors for the combined light curve. Table Q] presents the mid-transit 
times derived from fitting the models, as explained above, in the four observed 
transits for HAT-P-23, and Fig. [T] (bottom with diamond symbols) shows the 
resulting combined light curve for HAT-P-23. 



3.2. Ephemeris Determination 



We derived an improved orbital period for the HAT-P-23 system by per- 
form ing a least-squares linear fit to our data and the discovery ephemeris 
from iBakos et all ( 2011 ). weighing the individual m i d-tran sit times by their 
uncertainties. We have converted the IBakos et al.l (|2011l ) Barycentric Ju- 
lian Date based on Coordinated Universal Time (BJD_UTC) to the improved 
Dynamical Time-based system (BJD_TDB) as suggested, for example, in 
Eastman. Siverd fc Caudil t0l(h . The difference between UTC-bascd and 
TDB-bascd timings is a systematic offset which depends on recent additions 
of leap seconds to UTC. To convert BJDJJTC to BJD.TDB 0.00077 days 
must be added for transits observed after 2009 January 1st (JD 2,454,832.5). 
Our analysis yields an ephemeris for the system with a period of 1.2128868 ± 
0.0000004 days and an epoch of 2, 454, 852.26542 ± 0.00018 (BJD.TDB). The 
timing residuals are sho wn in Fig. [2] Our re sult improves the period reported 
in the discovery paper bv lBakos et al.l (|2011[ ) and also reduces the uncertainty. 
No transit-timing variations (TTVs) from a single period, which would suggest 
the presence of other planets in the system, can be claimed at this time. 



4. TRANSIT MODELING 



To estimate the essential transit parameters (the planet-to-star radius ra- 
tio: Rp/R*, the scaled semi-major axis: a/R+, and the orbital inclination 
of the system: i) from our combined ligh curve, we us e the Transit Anal - 
ysis Package (TAP) soft ware package des i gned for IDL (jGazak et al.1 120121) . 
It uses the now standard Mandel fc Agoll ( 2002 ) models for the transit light 
curve calculation and implements a Markov Chain Monte Carlo (MCMC) al- 
gorithm to find the best-fit parameters for the observed light curve. This 
analysis evaluates different combinations of parameters until it converges to 
an optimum solution. The packa ge also implements the wavelet-based noise 
treatment algorithm developed bv lCarter fc Winnl (|2009f ) . Such a scheme also 
allows estimating the uncertainty of the transit parameters. The Monte Carlo 
metho d implemen ted in TAP is a Metropolis-Hastings algorithm (see for ex- 
ample |Ford| ( 2006 )) that finds the set of parameters that minimizes the value 
of x 2 m the analyzed light curve. 
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The program starts with an initial estimation of the parameters and modi- 
fies a parameter of the set randomly, within reasonable limits of the parameter 
space. If the new set is a better estimation, the program accepts it. If not, the 
new set can be accepted or rejected depending on a condition that is evaluated 
randomly. Each set of parameters is defined as a link, and after execution, the 
algorithm returns a chain of all the sets of parameters that were evaluated. In 
order to estimate the best-fit parameters, TAP performs a Bayesian inference 
on the resulting chains, provided that a state of convergence has been reached 
in each chain. All the chains are added together and, for each parameter, TAP 
calculates the median, and the 15.9 and 84.1 percentile levels. The median 
is reported as the estimation of the parameter and the mentioned percentile 
levels are used to estimate the la value. 

The TAP software, however, does not perform a joint fit of the photometry 
data with the radial velocity data. This is a limitation in this case since the 
value of a/R* also depends on the radial velocity parameters when the orbit 
is non-circular. The intrinsic geometrical parameter that comes out of the 
transits alone is the duration of th e transit, equ i valen t to (reciprocal of 

the half duration of the transit) in lBakos et al. J 2011). T he relation between 
a/R* is expressed by equation (1) in iBakos et al. ( 2011 ). and from this it is 
evident that there is a subtle dependency of a/i?* on the eccentricity and 
orientation of the orbit. Therefore, one cannot derive a/i?* without knowing 
these orbital parameters. The dependency vanishes for circular orbits, but in- 
creases with eccentricity. This is the reason we adopted eccentri city e and the 
longitude of pcriastron uj values reported in IBakos et al. ( 2011 '). The orbital 
period, also an important parameter for determinig a/R+, can be determined 
from radial velocity curves as well as by successive observations of transits, as 
explained in § 13.21 

To start de MC MC algorithm, we set the orbital parameters of the sys- 
tem as reported by IBakos et "all toilb (e = 0.106 ± 0.044, u = 118° ± 25°) 
except for the period, where we used our result from th e previous se ction, 
and quadratic limb- darkening coefficients for the star from lcia'ret ( 2000h . For 
the case of HAT-P-23 ( T Rff = 5905 ± 50, [F e/H] = 0.15 ± 0.04, and logg = 
4.33 ± 0.06 — also from IBakos et al.l (|2011l) ) the coefficients are a = 0.324 
and b = 0.339 for our photometric bandpass. Southwortbl ( 20081 ) concluded 
that there is no significant difference in using either a single linear (u) or 
two quadratic (a & b) coefficients to describe the stellar limb-darkening in 
the analysis of high quality ground-based data. We ran 10 chains with 10 6 
samples in each chain. The results o f our modeling are presented in Table [2] 
along with the original parameters of IBakos et al.l (|201ll ) for comparison. 

To verify our results we have performed alternative model solutions. One 
alternative consisted of doing a simultaneous fit to the four individual light 
curves (instead of the combined one) searching for a global solution to our 
three main model parameters (R p /R+, a/i?*, and i). In this case we set the 
period and mid-transit times to the values obtained in § [3] and used the light 
curves where the baseline had been corrected using the airmass function as 
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outlined in § [2j Not surprisingly the results and uncertainties were very nearly 
the ones obtained originally given the random nature of the Monte Carlo fits. 
This, however, outlines the fact that our uncertainties may be optimistic in 
the sense that the uncertainties derived from calculating the mid-transit times, 
the period, and the airmass correction to the baseline are not carried over to 
the Monte Carlo solution. In order to obtain more realistic values to our final 
model uncertainties we need to solve simultaneously for these parameters as 
well. In this second analysis we again use the TAP software to search for 
a global solution to our three main model parameters, and now also let the 
software find the best values for the mid-transit times of the individual light 
curves (T c i, T C 2, T C 3, T C 4), a baseline corrections for each one (though we 
are limited in the sense that TAP uses a linear fit and not one based on 
the airmass value of the observations), and the system period. After longer 
computational time the resulting model values are similar to the ones obtained 
in our original analysis, but the formal uncertainties have nearly doubled 
for the model parameters and tripled for the mid-transit times. These are 
included in Tables [1] and [2] for comparison. 

The differences in the uncertainties between the two methods are expected. 
A "true" global fit would require all model parameters to be solved simul- 
taneously with as much data as possible. This would also have to include 
additional parameters like the limb-darkening coefficients and the rest of the 
orbital elements of the system. Since our available data is limited we elected 
to fix as many of the parameters as possible using various methods before 
solving for the three parameters of interest. Our period determination, for 
example, is more accurate than the one that could be derived from our four 
light curves alone since we use a longer time baseline. By the same token, 
our mid-transit times are also optimized to the individual light curves and 
solved for using robust methods. Individual fitting of the light curves may 
also be warranted since simultaneously solving for a single period may not 
describe the system accurately in the presence of possible mid-transit tim- 
ing variations due to the presence of another planet in the system. We also 
feel that our airmass-dependent light curve baseline correction is a stronger 
model than a simple linear fit. Our guiding principle was the optimization of 
the individual parameter fits by minimizing the resulting uncertainties. After 
constraining those parameters as best as possible using independent methods, 
then we proceeded to solve for the remaining ones. 



5. DISCUSSION 



The results presen ted in Table [2] agree in general with those found initially 
bv lBakos et al.l (|201ll) . In some instances, we managed to decrease the uncer- 
tainty levels. In particular, our modeled inclination (i) and scale d semi-major 



axis ( a/R*) for the system agrees with the results obtained bv iBakos et al 



(|201ll ) within the la level of the measurement s. While our inclin ation uncer- 
tainty could be considered slightly larger than IBakos et al.l ([201 if ) , our scaled 
semi-major axis uncertainty is reduced by at least a factor of two. We obtain 
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a more central transit and a slightly larger distance of the planet from the 
star. In general, there is degeneracy in these two variables since a smaller 
impact parameter (an inclination closer to 90°) results in a longer transit du- 
ration, which is also the result of a larger semi-major axis. This effect can 
also be mimicked by a planet further away from its star (thus going slower) 
and having a larger impact parameter (crossing less surface area of the star) . 
Limb-darkening effects observed in light curves taken at different wavelengths 
could help us sort out which case better represents this system. This comes 
about because stellar limb-darkening in general decreases with the inverse of 
the wavelength, yielding "flatter" depths in the light curves and "steeper" 
ingress and egress profiles. Our results do not show the sense of this degener- 
acy which leads us to believe that they may be real, although observations of 
transits at other wavelengths, preferably in the mid-infrared (J-, H-, K-band) 
are still desirable to settle this question. The major discrepancy in our results 
is the scaled planet-to-star radiu s ratio (Rp/R*) w hich is smaller in our case 



and with similar uncertainties as iBakos et al.l (| 2 1 lh . This is evidenced by a 
less deep transit depth in the light curve. Given a stellar radius, our planet is 
~ 5.5% smaller, but wc must point out that our results would still overlap at 



the 3ct level of uncertainty. This agrees with the expectation of iFortnev et al 
(12008) which p r edicte d a smaller planet size (~ 8.4% smaller) than observed 



bv lBakos et al.l (|201l[ ) from their theoretical models for a planet of this mass. 



6. CONCLUSION 

We have shown in the present work that several exoplanet transit light 
curves obtained with small telescopes can be successfully combined to produce 
a higher quality transit light curve that can then be modeled and compared 
with observations performed with telesc opes of larger apert ure. Our results in 
general compare favorably with those of IBakos et al.l (|201ll ). with instances of 
improved uncertainties. In particular, we observe a smaller planet size which 
agrees better with the expectations of lFortnev et al.l (|2008l ). although it must 
be emphasized that there are still large uncertainties in planetary models and 
no firm conclusion can be obtained from the results presented here. However, 
we could not improve on the uncertainty regarding the inclination of the sys- 
tem, which would have helped answer the question regarding the projected 
angle between the o rbital plane and the stellar equatorial plane outlined in 



Moutou et al.l (|201ll ). Furthermore, the individual light curve observations 
are valuable in the sense that their mid-transit times can be accurately de- 
termined and this in turn would help refine the orbital period of the system 
and/or determine possible transit timing variations that could hint at the pres- 
ence of another exoplanet in the system. An accurate orbital period is also 
valuable because it can be compared with the timings of secondary eclipses 
of the system. Comparing these values gives a more accurate measure of the 
eccentricity and argument of perihelion of the system than the radial velocity 
method. HAT-P-23b in particular is an inflated hot Jupiter-type planet in 
a close-in eccentric orbit that would show relatively deep secondary eclipse 
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light curves that may be observed with ground-based instruments, and thus 
an accurate period is desirable for this system. Further observations of tran- 
sits and possible secondary eclipses at various wavelengths with either larger 
telescopes or combined smaller telescopes are desirable in order to improve 
on the model parameters and better characterize this important inflated hot 
Jupiter system. 

The authors are grateful to an anonymous referee for the insightful com- 
ments on this work. 
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Fig. 1. Individual light curves (^Sp fESff wftte a alferisks) observed through an Rc- 
band filter at UDEM Observatory, and final combined linghtcurve (bottom one with 
diamonds) with 2-minute bins. Superimposed is the best-fit model light curve used 
to determine the mid-transit times. 



FABLE 1 

MID-TRANSIT TIMES FOR HAT-P-23 




BM3 


Global Fit 


UT Date 


Mid- Transit Time 


Mid- Transit Time 




BJD.TDB a 


BJD.TDB a 


3 June 2011 

4 August 2011 
16 August 2011 
21 August 2011 


5715.84150 ± 0.00057 
5777.69845 ±0.00055 
5789.82574 ±0.00051 
5794.67860 ±0.00051 


5715.84202i°;E 

5777.698001^0153 
5789.826121™°^ 

5794.67746l°;K 



a Barycentric Julian Date based on Dynamical Time 
(2,450,000+). To convert to BJD.UTC (Coordinated 
Universal Time) subtract 0.00077 days. 
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Fig. 2. Timing residuals for HAT-P-23. Dia monds represent o ur observations 
and the square is the original ephemeris from fB akos et al.l l201lh . The weighed 
least-squares fit yields a period of 1.2128868 ± 0.0000004 days and epoch of 
2, 454, 852.26542 ± 0.00018 (JDB.TDB). 



TABLE 2 

MODEL SYSTEM PARAMETERS FOR HAT-P-23 



Parameter 


Combined 
Light Curve 


Global Fit 


Bakos et al. (2011) 


i 




87.2t^ 


85.1 ±1.5 


a/R* 


a OO+0.06 a 
^• zo -0.12 


A QO+0-11 a 
^■ oz -0.17 


4.14 ±0.23 


R P I R* 


u.iiuo_ 0013 


u - iu0 ' -0.0026 


0.1169 ±0.0012 



a This parameter was obtained using e and u> values from lBakos et al.l (|2011 



